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Abstract 

hi this talk, we will demonstrate some concepts of a simultaneous treatment of neutrino 
decays and neutrino oscillations in an illustrative manner. This includes topics such as 
phase coherence discussions and time delay effects of massive supernova neutrinos. 

1 Introduction 

Neutrino decay in vacuum has often been considered as an alternative to neutrino oscillations 
{e.g., in ||l|-^). Either neutrino decay only (especially for atmospheric or solar neutrinos) or 
sequential combinations of neutrino oscillations and decays (especially neutrino oscillations in 
matter followed by neutrino decay in vacuum: MSW-mediated/MSW-catalyzed solar neutrino 
decay) have been studied. However, simultaneous neutrino decays and oscillations are also a 



possible scenario |T0|. It involves several quantum field theoretical issues such as phase coher- 
ence |]TU|,|T^ . In this talk, we will show some peculiarities coming from this sort of discussions, 
introduced in a quite conceptual manner and illustrated by several examples. 



2 Majoron decay as an example 

In order to demonstrate several kinematics and coherence issues of a decay model, we choose 



Majoron decay as an example |jT2|-|T^. Let us assume a generic effective interaction Lagrangian 
such as 



i j 

where J is the Majoron field, Ui are Majorona mass eigenstates, and Qij are the Majoron coupling 
constants. First of all, we observe that only mass eigenstates and not fiavor eigenstates may 
decay. Second, decay into active as weh as sterile neutrinos is, in principle, possible with 
this type of Lagrangian. Third, we assume the secondary active neutrinos to be, in principle, 
observable, but the Majorons, to first order, not. 
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2.1 Re-direction of neutrinos by decay 



Re-direction of neutrinos by decay is a purely kinematical effect. Since (at least) a third particle 
is involved in a decay process, the secondary neutrino may slightly change direction because 
of energy and momentum conservation. Let us now investigate the consequences for neutrino 
beams and radially symmetric point sources. 



Neutrino beams 

Figure shows the geometry of a neutrino beam produced at S and directed towards the 
detector D with an intermediate decay at X. One can derive from the kinematics of Majoron 
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Figure 1: A neutrino beam produced by the source S and detected at the detector D, where 
D is also the area of detection. The decay may happen at X, changing the direction of the 
secondary neutrino by an angle 9. 

decay that the angle 6 is limited by a maximum angle ^max for z/j — > Uj decay 

^ - ^ Ami = - > 0. (2) 



The angle ^max is determined by the ^-dependence for a not too hierarchical mass spectrum. 
Thus, we obtain for relativistic neutrinos ^max ^ 1- One can show that active secondary neutri- 
nos are, in principle, observable for accelerator, atmospheric, and reactor neutrinos (Amf^ ~ 
assumed) [ffO . 



Radially symmetric neutrino sources 

From Fig. ^ we see that for decay of one parent neutrino into exactly one secondary neutrino 
the overall flux of a radially symmetric neutrino source is conserved. Thus, active secondary 
neutrinos are, in principle, observable for solar and supernova neutrinos. Furthermore, we 
notice that especially for massive supernova neutrinos, the travel times on different paths may 
be different even for small angles ^max- As we will see later, this will modify the time dependence 
of the signal at the detector. 



2.2 Interference effects 

In this section, we will study phase coherence in decay processes. Especially, we are interested 
in the observability of interference effects with intermediate decays between production and 
detection. Neutrino oscillations of decay products will be one example for such an interference 
effect. 
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Figure 2: A radially symmetric neutrino source and several paths of neutrinos with one 
intermediate decay. Drawing any concentric sphere around the source illustrates that the 
overall flux through any of these spheres is conserved as well as radially symmetric. 

Neglecting the Majoron field as well as the operators, we know that to first order in the 
S-matrix expansion 

2 

' Ant ~ ^fl'ij'^i'^i- (3) 

Thus, the interaction destroys an in state and creates an out state by application of the appro- 
priate annihilation and creation operators in the field expansions within the Lagrangian. Let 
us now assume an incoming and outgoing superposition of mass eigenstates, i.e., active fiavor 
eigenstates (ignoring neutrino propagation): 

R = Wa) = J2^:M), (4) 

i 

\out) = \u(,) = J2u;^\u,). (5) 
j 

Applying this to Eq. (^) shows that the differential decay rate may indeed contain interference 
terms oc {i^jlSijlui)* {uilSkili'k) with {i,j) ^ {k,l), corresponding to interference of different 
decay channels. 

From a different point of view, interference of different decay channels corresponds to co- 
herent summation of amplitudes. This is illustrated in Fig. ^ for the case of Majoron decay. 
Let us compare this to the examples of weak interaction processes shown in Fig. ^. In the case 
of W decay, the coherence among the out neutrinos of different flavors is destroyed by the mass 
differences of the participating leptons corresponding to the different flavors. In other words, 
the wave packets of the different neutrino flavors do not sufficiently overlap because of the 
kinematics of the different leptons. Thus, the Feynman diagrams for different flavors need to 
be summed incoherently. However, in the case of Z° decay, the produced neutrino-antineutrino 
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Figure 3: Incoherent (left) and coherent (right) summation of amphtudes for the case of 
Majoron decay. 
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Figure 4: Incoherent (left) and coherent (right) summation of amplitudes for some cases of W 
(left) and Z° (right) decay. 

pairs have small enough mass differences to allow wave packet overlaps of different flavor or 
mass eigenstates. This means that we need to sum the Feynman diagrams coherently [[l^] . 

Comparing decay to Majoron decay, we conclude that interference is not necessarily 
destroyed in neutrino decay. Thus, depending on the decay scenario, we may have to take into 
account neutrino oscillations of secondary neutrinos or parent neutrinos, such as illustrated in 
Fig.i 




Figure 5: Examples of decay scenarios in which neutrino oscillations of secondary neutrinos 
(left) or neutrino oscillations of parent neutrinos (right) become, in principle, possible. For 
these scenarios the Majoron coupling constants shown in the figure need to have non-zero 
values. The circles correspond to the mass eigenstates z/3, 1/2, and z/i, respectively. 



2.3 Decay as measurement 



Similar to the production and detection processes in neutrino oscillations, neutrino decay acts 
as a measurement. In general, decay destroys an incoming superposition of mass eigenstates 
and creates an outgoing one, as it is indicated in Fig. |[ If we can only detect the secondary 




Figure 6: Decay process illustrated, with an in and out superposition of mass eigenstates. 

neutrinos but not the Majorons, we will in many cases (for unchanged quantum numbers and 
similar energies) not even be able to tell, if there has been a decay between production and 
detection, or not. However, since there is a third particle involved (the Majoron), these two 
cases can, in principle, be distinguished. This is equivalent to the fact that Feynman diagrams 
of different orders do not interfere. Thus, decay acts as a measurement and resets the relative 
phase among the mass eigenstates in the out state to 0, similar to the example of decay 
above. Therefore, for any considered neutrino oscillation (before or after decay) the oscillation 
phases at the detector depend on the position of the decay. Figure |^ shows the number of 
neutrinos over the traveling distance for small and large decay rates. For small decay rates. 
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Figure 7: Fraction of undecayed neutrinos N/Nq over the fraction of the baseline l/L for 
small (left) and large (right) decay rates (neutrinos treated as particles). For exponential decay 
oc exp{—al/E) this corresponds to a <^ E/L (left) and a ~ E/L (right). 

the positions of decay are almost equally spread over the entire traveling distance. Thus, any 
oscillation phase will be averaged out, similar to the case of production or detection regions 
larger than the oscillation length. However, for large decay rates more neutrinos will decay in 
the beginning of the path than at the end. Since the oscillation phases are averaged over all 
possible decay positions, we may thus expect a net oscillatory effect. 
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3 Invisible decay products 



For decay into unobservable particles, such as sterile decoupled neutrinos, we do not have to 
take care of secondary neutrinos as well as the type of neutrino source. Therefore, this is the 



simplest case of neutrino decay. One can show that the transition probability is given by |T0 



with 



Pcp conserving 

2 ^53J^^sin2Ai,e-^'^ (6) 



i>j 



Jff = U.^U:^U;,U^„ A,, = ^^, (7) 

Here ctj = mi/ro^i, where ro,i is the (rest frame) lifetime of Ui. In Eq. (|^) we see that the 
oscillatory terms are damped by exponentials describing the disappearance of neutrinos into 
decay products invisible to the detector. 



Example: Atmospheric neutrino decay 

In order to demonstrate the effects of invisible decay, we may choose a decay scenario similar 
to the one in Ref. |0 shown in Fig. | for atmospheric neutrino decay. For the parameters we 




Figure 8: A scenario for atmospheric neutrino decay. Here the mass eigenstate z/2 may decay 
into a sterile eigenstate, which is not mixing with the active mass eigenstates. 

take cos^6'23 = 0.30 as well as Amgg = 3.3 ■ 10^'^ eV^, since we want to take into account 
neutrino oscillations in addition to neutrino decay. Figure || shows the survival probability P^^ 
for different decay rates a and sensitivities L/E. One can imagine the transition curve between 
decay and oscillation dominated regions, determined by aL/E ~ 1. The purple line in Fig. ^ 
corresponds to the cut shown in Fig. ITUI, in which the damping of the oscillation by decay can 
be clearly seen. 
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Contour plot of the survival probability P^^ over the decay rate a and the sensitivity 
Top. The decay scenario used is shown in Fig. ^ For the 
atmospheric neutrino oscillation parameters we take cos^ 623 = 0.30 and = 3.3 ■ 10"'^ eV^. 



Figure 9: 
L/E (similar plot as in Ref. 
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Figure 10: Cut through Fig. ^ at a = 1/1000 GeV/km at the purple line. 

4 Visible decay products in neutrino beams 

We showed for the example of Majoron decay that decay into active neutrinos, or sterile neu- 
trinos mixing with active ones, involves more complicated discussions than decay into invisible 
neutrinos, such as sterile decoupled neutrinos. In this section, we will only give a notion of the 
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results for visible secondary neutrinos in neutrino beams. 

Let us define P^^ to be the transition probabilities for the flavor transition — >■ 1^/3 with 
exactly i intermediate decays {i = 0, 1, 2, . . .). In order to calculate the total transition proba- 
bility, the transition probabilities for different indices i have to be summed over or not. This 
depends on the ability to distinguish the secondary from the parent neutrinos, i.e., conceptual 
properties of the detector and the problem. For example, for decay into antiparticles the decay 
products and the parent neutrinos may have different signatures in the detector. Another ex- 
ample is energy resolution: since neutrinos loose some energy to third particles by decay, the 
detector may distinguish the parent and secondary neutrinos by its energy resolution. 

Assuming that no secondary neutrinos escape detection by kinematics, one can show for the 
first transition terms P^^ that [lU 



pO 



= P 



invisible 

a/3 5 



(9) 



EE 

ij kl 



X |3?(ir:^.^J [(r,, - r,,) (e-r- cos(2Afc,) - e-^^' cos(2A,,)) 

- 2 (Ai,- + Aifc) (e-^''= sin(2Afe,) - e"^^' sin(2Ai,-))] 

- [(r,, - r,fc) (e-r- sin(2Afc,) - e~^^^ sin(2Az,)) 

+ 2 (Ai,- + A;fc) (e-^"= cos(2Afci) - e"^^' cos(2Ai,))] | 



(10) 



p2 

-'a/3 



where P^^^'^^ibie 



is given by Eq. (H), -ft'"^; = U*^Uf3jUakU^i is a generalization of J^"^, and atj 



f^i I To^ij is the decay rate for the channel z/j 



analogously defined to at = rrii/To^i = 



Example 

In order to show the effects for visible decay products, we construct an example with the decay 



scenario shown in Fig. 11. Let us look at the survival probability of electron neutrinos Pee and 




Figure 11: Example for a decay scenario with > 1712 > mi. 



ignore higher order decay effects, i.e., only consider Pg° and Pgg. For the decay scenario in 
Fig. |Il] one can split up the non- vanishing terms in the sum in P^^ in Eq. (|T0|) into 



P^ 
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where P^^ describes the production of new by decay and -Peg'™* the interference effects 
(neutrino oscillations) before or after decay. Figure shows the separated signals and their 
sum for the parameters in the figure caption. The terms P}^^-, describing the production of new 
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Figure 12: The separated survival probabilities Pl^ for decay into individual mass eigenstates 
z/fc, the interference probability P}^^^-, as well as the total survival probability with exactly one 
intermediate decay P}^ (similar plot as in Ref. ||10[)- Here the decay scenario in Fig. ^is used 
with parameter values Am|2 = 3.3 ■ IQ-^eV^ Am|i = 5 ■ 10-^eV^ ^2 = 1/1000 GeV/km, 
03 = 1/1100 GeV/km, as well as trimaximal mixing. 



mass eigenstates by decay, are exponentially growing in the beginning. For large L/E, P^f is 
falling again because i>2 decays with a larger rate than being produced. The interference term 
p^i,int jg basically determined by the two beat frequencies induced by the two Am^ involved, 
i.e., neutrino oscillations before and after decay. 

Taking P^^ into account, which describes the survival probability of the neutrinos arriving at 



the detector without any decay between production and detection, yields the result in Fig. 13 



Note that and P^, may only be sensibly added, if the detector cannot distinguish between 
parent and secondary neutrinos. 



5 Decays of supernova neutrinos 

In this section, we will focuse on supernova neutrinos. We will show certain properties of 
supernova neutrino propagation as well as the effects modifying neutrino event rates. 



5.1 Issues especially concerning supernova neutrinos 

Since a supernova may be approximated as a far-distant point source, we have to incorporate 
some new concepts {e.g., Refs. ||lT],[l^): 
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Figure 13: The survival probabilities without intermediate decays P^^ and with one inter- 
mediate decay P}^ between production and detection, as well as their sum (similar plot as in 
Ref. [lO]). The parameter values are chosen as given in the caption of Fig. O. 



Time delays of massive neutrinos 

For massive neutrinos the velocity of propagation depends on the mass. Even on the direct 
path L, massive neutrinos will be delayed by At ~ Lm? / {2E'^). In addition, re-direction by 
decay opens the possibility for paths from the supernova to the detector other than the direct 
path, such as shown in Fig. |l^. Again, neutrinos will be delayed by an additional time interval. 




Figure 14: A traveling path from the supernova SN to the detector D different from the direct 
path along the baseline L. 

though the change of direction by decay is in many cases quite small. 

In order to investigate this effect, let us assume a radially symmetric source producing 
neutrinos of flavor z/q, at t = — L, i.e., 

K\t) = '^-^ = N^6{t + L\ (12) 

so that 

oo 

j ^':\t)dt = N^. (13) 

— oo 

Therefore, for such a source flux massless neutrinos would arrive at t = 0. Let us further define 
$^A*(t) to be the number of neutrinos per time interval at the detector, which are produced 
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as flavor z/^ and detected as flavor with exactly i intermediate decays. This deflnition is 
completely analogous to P^^, but in addition takes into account the time dependence of the 
signal. 

Loss of coherence because of long baselines 

We know from the wave packet treatment of neutrino oscillations that the coherence length of 
neutrino oscillations with Am^^ is given by [p!^-pT| 



- coh 



(14) 



Here 



a: 



J\2 



(15) 



is the combined wave packet width a^. of the production P and detection D processes. The 
decay process X can be treated similarly to an intermediate process between production P and 
detection D by using 



+ or iaif 



(„X\2 I („D\2 



(16) 



depending on what processes we are looking at. 



5.2 Dispersion by different neutrino masses 

As a flrst effect, which we illustrate with an example, we demonstrate the propagation of mass 
eigenstates traveling with different velocities before and after decay because of the different 
masses of parent and secondary neutrinos. Thus, the arrival time depends on the position of 



decay. We choose the decay scenario in Fig. |I5| and assume incoherent propagation at all times, 
i.e., the travel distance between any two processes in the problem is much longer than the 
respective coherence length. In addition, we ignore effects of different traveling path lengths as 




Figure 15: A decay scenario with ms > m2 > mi for decay of one mass eigenstate into two 
different decay channels. 

well as repeated decays. Figure |16| shows the results for the parameter values given in the flgure 
caption. For no intermediate decays between production and detection, the source pulses are 
also detected as pulses. For one intermediate decay, the neutrinos travel with the velocity of the 
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Figure 16: The signals for no intermediate decays between production and detection 

as well as for one intermediate decay separated into the terms for ui and z/2 as decay 



products (similar plot as in Ref. [11 
parameter values 
"32 = "31 = E/L, 



In this example, the decay scenario in Fig. 15 with the 
4eV, ma = 2eV, mi = leV, L = lO^^m ~ 32kpc, E = lOMeV, 
9 ■ 10^ AttL'^/D, as well as trimaximal mixing is used. Note that in this 



example the flavor index a may refer to any flavor. 



heavy parent mass eigenstate before decay and with the velocity of the light secondary mass 
eigenstate after decay. Thus, one can find the exponential distribution of the decay positions 
in the time structure of the signal at the detector. This is an effect which may affect or even 
mimic the signal structure expected from supernova models. 

5.3 Dispersion by different traveling path lengths 

In this section, we will illustrate the effects of different traveling path lengths with an example. 
Note that the smaller the maximum re-direction by decay 6'niax is, the smaller the dispersion by 
time delays on different traveling paths becomes (c/.. Fig. |T^. However, there will be dispersion 




Figure 17: Different paths from the supernova SN to the detector D, as well as the path with 
maximum length by kinematics and geometry (red). 



by different neutrino masses even for small 6'max, such as it was shown in the example in Sec. 



We use an example with the decay scenario shown in Fig. 18 and the parameter values from the 



last example in Sec. ^]2| (c/., caption of Fig. 0). In addition, we approximate the differential 
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Figure 18: Example for a decay scenario with > m2 > mi. 



decay rate by its mean 



dr 



dcosO 



1-COS Sniax 





for 

otherwise 



[17) 



with an effective 6'max < ^max- Thus, we may choose ^max smaller than the actual 6*^0x5 in order 
to investigate the dependence on the path lengths. Figure |19| shows the (approximated) signal 
with one intermediate decay for several values of ^max- The figure demonstrates that for large 
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Figure 19: Approximated signal with one intermediate decay for several values of 6'inax defined 
in the text (similar plot as in Ref. [11 1). For the decay scenario we use the one in Fig. |l^. The 



parameter values are chosen such as in the last example in the caption of Fig. ITB 



6'inax late time arrivals are favored and early time arrivals suppressed. For small 6^ 
observe an exponential behavior such as expected from the last example. 



we almost 



5.4 Early coherent decays 

For this effect, we assume the decay rates to be large enough such that the neutrinos are 
still coherently propagating at the position of decay, but loose coherence between decay and 
detection. This also implies that all neutrinos decay before detection. 

For the decay scenario we need to have simultaneous coupling of two mass eigenstates to 
the decay product, such as in Fig. [2^. In addition, we assume trimaximal mixing. Since it can 
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Figure 20: Example for a decay scenario with ms > 1712 > mi and two decay channels into the 
same decay product. 

be shown that the detector can in most cases not resolve the time dependence of the signal, we 
integrate the flux over time: 

00 

A^3' = / ^a,'dt. (18) 
—00 

Moreover, similar to Eq. (|TT]), we define N[^t to be the number of neutrinos coming from 
interference terms in the calculation and A^mcoh to be the number of neutrinos coming from 
incoherent propagation. Figure |2T] illustrates that for small Am^j <^ a interference effects 
become most important. For larger Am|2 more and more oscillations take place until the 
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Figure 21: The ratio of neutrinos coming from the calculation of interference terms A^int to 
the number of neutrinos coming from incoherent propagation A^incoh- This ratio is plotted over 
the ratio of Am32, which causes the interference terms, to the decay rate a. For the calculation 



trimaximal mixing as well as the decay scenario in Fig. 20 is assumed 



neutrinos decay and thus the oscillation phases at the position of decay are more and more 
averaged out over all possible decay positions. Note that we may observe an interference effect, 
even if only one stable neutrino arrives at the detector. 
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6 Summary 



So far, parameters have only been fitted for neutrino oscillations or special neutrino decay 
scenarios. However, fitting all oscillation and decay parameters may produce new solutions. 
Furthermore, decay effects for supernova neutrinos have been ignored. Taking them into ac- 
count may alter or even mimic the signals expected from supernova models. Finally, supernova 
neutrino observations have only indicated that at least one mass eigenstate is stable. Neverthe- 
less, interference phenomena may have to be taken into account, even if only one stable mass 
eigenstate arrives at the detector. Since non-zero values for neutrino masses imply, in principle, 
not only neutrino oscillations, but also neutrino decay, we conclude that neutrino decay should 
be incorporated into the general neutrino oscillation discussion. 
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